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■ Theory predicts and observations confirm that low-mass stars (like the Sun) in their 



early life grow by accreting gas from the surrounding material. But for stars ~ 10 times 
more massive than the Sun (~ 10 Mq), the powerful stellar radiation is expected to 
inhibit accretion^ and thus limit the growth of their mass. Clearly, stars with masses 
> 10 Mq exist, so there must be a way for them to form. The problem may be solved by 
non-spherical accretion^ '^, which allows some of the stellar photons to escape along the 
symmetry axis where the density is lower. The recent detection of rotating disks^^^ and 
toroids"^ around very young massive stars has lent support to the idea that high-mass 
8 Mq) stars could form in this way. Here we report observations of an ammonia line 
towards a high-mass star forming region. We conclude from the data that the gas is 
falling inwards towards a very young star of ~ 20 Mq, in line with theoretical predictions 
of non-spherical accretion. 
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In the work we report here, we have revealed the simultaneous presence of three elements in 
the same massive object: outflow, rotation, and infall. Although evidence of circumstellar disks or 
toroids in massive objects had been provided by a number of observational studies*^" (see Ref. 11 
for a review), it remained to be proved that these are accretion structures, sustaining the growth 
of the central star(s). Evidence of infall had been found only in a very limited number of massive 
objects^^"^^, but in no case had the simultaneous presence of rotation, outflow, and infall towards an 
0-B (proto)star been established. Our findings hence represent a substantial advance in this field. 

The source under investigation is G24.78+0.08, a massive star forming region located at a distance 
of 7.7 kpc, studied by us in great detaiF'^^"^^. This contains a cluster of Young Stellar Objects 
(YSOs), three of them associated with rotating toroids. In one of these, G24 Al, the presence of an 
early-type star is witnessed by a hypercompact Hll region (M.T.B. et al, manuscript in preparation) 
with a diameter ^Cf.'2 (^1,500 AU) and located at the geometrical centre of the toroid (see Fig. 1). 
The continuum spectrum resembles that of a classical (Stromgren) Hll region around a zero-age main 
sequence star of spectral type 09.5^^, corresponding to a mass of ~ 20 Mq, a luminosity of ~33,000 Lq 
and a Lyman continuum of 5 x 10^^ sec~^. Being the Hll region very bright at centimeter wavelengths 
(^2,000 K), it is easy to observe the colder molecular gas (--^100 K) in absorption against it. If the 
toroid is not only rotating, but also accreting onto the central star, one expects to see absorption at 
positive velocities (i.e. red-shifted by Doppler efi^ect) relative to the stellar velocity The latter 
(~110.8 km s~^) is equal to the mean velocity of the gas in the toroid, estimated from molecules whose 
lines are not affected by absorption^^, such as methyl cyanide (CH3CN). 

With this in mind, we observed simultaneously the continuum emission and the ammonia (NH3) 
(2,2) inversion transition at 1.3-cm wavelength using the Very Large Array (VLA) interferometer of 
the National Radio Astronomy Observatory (NRAO) in the B configuration. As expected, the line 
is seen in absorption towards the continuum of the Hll region: this is shown by the white contours 
in Fig. la, where one can also see the emission (black contours) detected towards the adjacent toroid 
G24 A2, not discussed here. Clearly, the peak of absorption as well as the peak of the 1.3-cm continuum 
emission (i.e. of the Hll region) lie at the center of the toroid previously imaged by us (yellow line) and 
along the axis of the associated bipolar outflow (black arrows). The velocity gradient in the rotating 
toroid is shown in Fig. lb, overlaied on an image and contour map of the CH3CN emission at 1.4-mm 
wavelength previously obtained by us^'^^ with the IRAM Plateau de Bure interferometer (PdBI) in 
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the most extended configuration. 

Figure 2 illustrates the most important finding of our study. In order to understand this figure, 
one must keep in mind that the NH3 inversion transitions have a complex structure, made of a "main" 
component and four "hyperfine satellites". The former line is intrinsically stronger than the latter 
ones by a factor ~14, so that, roughly speaking, emission in the main line is representative of lower gas 
densities, whereas the satellites trace the densest material. In Fig. 2 we show position-velocity plots 
for both the main line (Fig. 2a) and the mean satellite emission (Fig. 2b), where the offset in position is 
measured along the plane of the rotating toroid. In Fig. 2c we show the same plot for the CH3CN(12- 
11) = 3 line which was used by us to trace rotation in the toroid^. Three considerations are in 
order: 1) Towards the hypercompact Hll region, the satellite absorption is strongly biased towards 
positive velocities, i.e. red-shifted with respect to the velocity of the star; 2) The NH3 main line is 
seen in absorption at both positive and negative velocities, but the blue-shifted absorption is fainter, 
broader, and more optically thin than the red-shifted one; 3) The velocity gradient seen in the CH3CN 
line (outlined by the green line in Fig. 2c), is detected also in the NH3 main line (green line across the 
two emission peaks in Fig. 2a), thus confirming the presence of rotation. 

The important conclusion is that the densest gas seen in the satellite absorption is moving away 
from the observer at a speed of <^2 km s"^ (see Fig. 2b) towards the hypercompact Hll region, namely 
towards the star. This fact proves the presence of infall onto the 0-type star at the center of the 
toroid. On the other hand, the faint and broad blue-shifted absorption seen in the main line is very 
likely due to the low-column density gas in the lobe of the molecular outfiow that is moving towards 
the observer: this is indicated by the lower main line/satellite ratio (i.e. lower optical depth) observed 
in the blue line wing with respect to the one in the red-shifted absorption. Incidentally, we note that 
the broad blue-shifted absorption can be used to discriminate between sources G24 Al and G24 A2 as 
origin of the outfiow. Since the blue-shifted lobe is directed towards NW^^, it cannot originate from 
G24 A2 or it would not pass in front of G24 Al. We hence believe that the star in G24 Al is powering 
the flow. 

What are the implications of the detection of infall in G24 Al? Using the properties of the 
NH3(2,2) transition and an estimate of the gas temperature^^, one can derive the column density of 
the infalling gas, A^Hj, from the ratio between the main line intensity and that of the satellites'^, 
assummg an NH3 abundance relative to H2 of 10-^-10"^ (refs 18, 23). From A^Ha and the infall 
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velocity, Vi^f — 2 km s , assuming free-fall one obtains the mass accretion rate onto the star inside a 
solid angle fi: Mace — r2/(47r) (4 x 10~^-10~^) Mqyt^^. A detailed derivation of the mass accretion 
rate is presented in Supplementary Information. The range of values reflects the uncertainty on N^^ 
and on the radius at which Vit^i is measured. Noticeably, Mace is much larger than the critical rate 
above which formation of an Hll region is inhibited^^ if accretion is spherically symmetric (i.e. for 
O = 47r). For an 09.5 star this is Mmh — 8 x 10~^ Mq yr"-*^, much less than Maee^ this should suffice 
to prevent the formation of an Hll region. The fact that, instead, an Hll region is detected can be 
explained only if the accretion is not spherically symmetric (O < 47r), which in turn supports the 
existence of a circumstellar disk-like structure, detected by us on a larger scale as a toroid. 

It is also interesting to estimate the rate needed for the ram pressure to overcome the thermal 
pressure of the Hll region (Mhii $7/(4 tt) 2 x 10~^ M0yr~^). This has been obtained assuming an 
electron density and diameter of the Hii region respectively of ~ 10® cm~^ ^^^^ and 0'.'2 (see above) or 
0.0075 pc. One sees that Mace > Mhii independently of 0,, which indicates that the accretion is large 
enough to brake or slow down the expansion of the Hll region through the toroid. This conclusion is 
strengthened if one takes into account the effect of stellar gravity, likely non-negligible in such a small 
Hll region. Detailed models support this conclusion, and predict that hypercompact Hll regions could 
be long lived because of trapping of the ionized gas^^ owing to infall in the gravitational field of the 
star. 

We now consider the effect of our findings on current theories of massive star formation. Stars as 
massive as ~ 8 Mq are believed to reach the zero-age main sequence still accreting. At this stage, 
their powerful radiation pressure appears to be strong enough to halt the infalling material^'^®, which 
should inhibit further growth of the star beyond the limit of ~ 8 Mq. In order to solve this problem 
merging of low- mass stars has been proposed^'', although - even in the most favourable case of induced 
binary mergers - to be effective this scenario requires very large stellar densities, at least of the order 
of 10® stars pc~^. Increasing the mass accretion rate may be a viable theoretical solution to overcome 
the radiation pressure with the ram pressure of the infalling material. Moreover, this would solve the 
lifetime problem^^ related to the fact that with accretion rates typical of low-mass stars (10~^ Mq yr~^) 
the time needed to form a star > 10 Mq would be too long (> 10® yr), comparable to the entire lifetime 
of the star. Higher rates can be achieved either by assuming a large turbulent pressure^^'^^ in the 
molecular cores forming 0-B type stars, or deepening the gravitational well^° around a massive star 
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by surrounding it with a tight cluster of lower mass companions. Non-spherical accretion may also 
soften the problem^'^: focusing accretion into a circumstellar disk reduces significantly the radiation 
pressure force along the equator by beaming the radiation in the direction of the poles. 

In our case, direct accretion onto the central massive (proto-)star has not yet been detected, so 
that the system may be in a (relatively) late phase, when the final mass of the central object has 
been assembled and the infalling material is being deflected in the outflow instead of being accreted. 
Nevertheless, our detection of infall towards G24 Al is an important milestone towards an observational 
constraint of the theoretical models. In fact, even if the very accretion phase were finished, the presence 
of a rotating toroid infalling on a hypercompact Hll region, located at the base of a powerful collimated 
bipolar flow, is broadly in agreement with the expectations of the non-spherical accretion scenario for 
the formation of massive stars. Our results thus suggest that this is a plausible mechanism for the 
formation of stars at least as massive as 20 Mq. More 'exotic' mechanisms (such as stellar mergers) 
might still be required to form stars of even higher mass. 
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Figure 1: Absorption and emission by molecular gas towards the hypercompact Hll 
region G24 Al. a, Map of the intensity integrated under the NH3(2,2) inversion line (colour image). 

Black contours indicate positive intensity (emission), while white contours are negative (absorption). 
The white filled circle marks the position of the hypercompact Hll region detected in the continuum 
emission at 1.3-cm wavelength^^ (Beltran et al. in preparation). The black arrows outline the direction 
of the bipolar outflow, ejected along the axis of the disk. The disk plane is denoted by the yellow 
line. Note the deep absorption towards the Hll region, at the center of the toroid. The final maps 
have been reconstructed with a circular synthesized beam of O'.'S x O'.'S, which is represented by the 
hatched circle. The 1-cr noise of the integrated NH3(2,2) line is 0.6 mJybeam"^. Contour levels start 
at a ±3-cr level and vary by 3-a (positive contours) and 12-cr (negative contours), b. Map (image and 
contours) of the intensity integrated under the CH3CN(12-11) K = S line^^ outlining the toroids in 
G24 Al and G24 A2 (the latter is not discussed in this article). The colour scale corresponds to the 
velocity field of the toroid in G24 Al: note how the velocity changes gradually from red-shifted to 
blue-shifted, consistently with rotation about the outflow axis. A velocity gradient, similar to that in 
G24 Al, has also been detected^ in the other core G24 A2, but is neither shown nor discussed here. 
Symbols have the same meaning as above. The resulting synthesized beam is 1'.'2 x Of! 5 at PA = —174° 
and is represented by the hatched ellipse at the bottom right. The 1-a noise of the intensity integrated 
under the CH3GN(12-11) K = 3 line is 35 mJybeam"^. Contour levels start at 0.25 Jybeam"^ and 
increase by 0.2 Jybeam"^ . 6, declination; a, right ascension. 

Figure 2: Velocity field in the massive toroid G24 Al. a. Position- velocity plot of the 
NH3(2,2) main line intensity. The positional cut is made along the plane of the disk, PA=— 135°, 
namely along the yellow line in Fig. 1. Both offset and velocity are computed with respect to the 
star. Black and white contours correspond respectively to emission and absorption. The continuum 
has been subtracted from the line emission. The 1-a noise is 1 mJybeam"^. Contour levels start 
at a 3-(T level and increase by 3-a (positive contours) and 6-a (negative contours). The error bars 
denote the angular and spectral resolution. The original spectral resolution was ~ 0.31 km s~^, which 
was degraded by a factor 2 to improve the S/N by averaging pairs of adjacent channels. Note the 
velocity gradient outlined by the green line, b. As a but for the NH3(2,2) satellite line. Note how the 
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absorption peak is red-shifted with respect to the stellar velocity, c, As a but for the CH3CN(12-11) 
K = S line: the velocity gradient outlined by the green line is indicative of rotation^^ . The 1-cr noise 
is 48 mJybeam"^. 
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